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Abstract 

We report improved measurements of branching fractions for B — > Kir, 7r + 7r~, 7r + 7r° and KK 
decays based on a data sample of 85.0 million BB pairs collected at the T(45) resonance with the 
Belle detector at the KEKB e + e~ storage ring. This data sample is almost three times larger than 
the sample previously used. We observe clear signals for B — > Kir, it + it~ and ir + ir° decays and set 
upper limits on B — > KK decays. The results can be used to give model-dependent constraints on 
the CKM angle </>3, as well as limits on the hadronic uncertainty in the time-dependent analysis of 
the angle <p2- 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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Recent studies at B factories have significantly improved our knowledge of heavy-flavor 
physics. In particular, the establishment of mixing-induced CP violation in the 5-meson 
system |l|,|2j] is encouraging for further tests of the Standard Model based on determinations 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 

.B-meson decays to Kir, tttt and KK final states are dominated by b — > u tree and 
b — > s, d penguin diagrams. The properties of these decays provide information that can 
be used to determine the CKM angles 02 and 03 0]. However, the extraction of these 
angles suffers from hadronic uncertainties present in the current theoretical description and 
from the small amplitudes of b —>■ u, s, d transitions. To solve these difficulties, various 
theoretical approaches based on flavor symmetries and dynamical calculations in the heavy- 
quark limit p have been proposed. In order to utilize these methods, the precision of the 
existing experimental results p, 0, H, H, Ei EH must be improved. 

In this paper, we report updated measurements of the branching fractions for B — > 
Kit, 7T + 7r~, 7r + 7r° and KK decays. Recent results for B° — > 7r°7r° have been reported 
elsewhere |l3l Eij . The measurements reported here are based on a 78 fb _1 data sample 
collected at the T(45) resonance, with the Belle detector [l5| at the KEKB e + e~ storage 
ring [16]. This sample corresponds to 85.0 ± 0.5 million BB pairs and is about three times 
larger than that used for our previous analysis |6] . The previous results are superseded with 
significantly improved statistical precision. Throughout this paper, neutral and charged B 
mesons are assumed to be produced in equal amounts at the T(4S f ). The inclusion of the 
charge conjugate decay is implied, unless explicitly stated. 

The Belle detector is a large-solid-angle spectrometer consisting of a three-layer silicon 
vertex detector, a 50-layer central drift chamber (CDC), an array of threshold Cherenkov 
counters with silica aerogel radiators (ACC), time-of-flight scintillation counters, and an 
electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return located 
outside of the coil is instrumented to detect K\ mesons and to identify muons. A detailed 
description of the Belle detector can be found elsewhere (lif . 

The basic analysis procedure is the same as described in Ref. Q. However, the data sam- 
ple used in this analysis was reprocessed with an improved tracking algorithm that reduces 
the probability of incorrectly associating CDC hits in the track finding. This improvement 
changes the efficiencies for the kinematic reconstruction of the signal as well as for the mea- 
surement of specific ionization energy loss (dE/dx) in the CDC from the values given in 
Ref. @- 

The 7r ± mass is assigned to each charged track. Tracks used to form B candidates 
are required to originate from the interaction region based on their impact parameters. 
Kg mesons are reconstructed using pairs of oppositely charged tracks that have invariant 
masses in the range 480 MeV/c 2 < M n7T < 516 MeV/c 2 . A reconstructed Kg is required to 
have a displaced vertex and a flight direction consistent with that of a K s originating from 
the interaction region. Pairs of photons with invariant masses in the range 115 MeV/c 2 < 
M 77 < 152 MeV/c 2 are used to form 7r° mesons. The measured energy of each photon in 
the laboratory frame is required to be greater than 50 MeV in the barrel region, defined as 
32° < # 7 < 128°, and greater than 100 MeV in the end-cap regions, defined as 17° < # 7 < 32° 
or 128° <9 1 < 150°, where # 7 denotes the polar angle of the photon with respect to the e~ 
beam. Signal B candidates are required to satisfy 5.27 GeV/c 2 < M^ c < 5.29 GeV/c 2 and 



-0.3 GeV< AE < 0.5 GeV, where M bc = yfEg am -pg, AE = E* B — E* eam , £* eam is the 
beam-energy, and p* B and E* B are the momentum and energy of the reconstructed B meson, 
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all evaluated in the e + e~ center-of-mass (CM) frame. The signal efficiencies of the kinematic 
reconstruction, estimated using GEANT-based Monte Carlo (MC) simulations, are listed 
in Table |U 

Charged tracks from B candidates have momenta ranging from 1.5 up to 4.5 GeV/c in 
the laboratory frame. They are distinguished as or 7r ± mesons by the number of photo- 
electrons (iVp. e .) detected by the ACC and dE/dx measured in the CDC. These quantities 
are used to form a K ± identification (KID) likelihood ratio IZk = £k /(£%: + Ar)> where Ck 
denotes the product of the individual likelihoods of N p _ e _ and dE/dx for K mesons, and L v 
is the corresponding product for mesons. The requirements on IZk used in this analysis 
yield a identification efficiency of 84.4% with a 7r ± misidentification rate of 5.3% for 
K ± candidates, and a ir^ identification efficiency of 91.2% with a K ± misidentification rate 
of 10.2% for tt^ candidates. The efficiencies and misidentification rates are measured by 
comparing the yields of high-momentum D* + -tagged D° — ► K~n + decays before and after 
applying the IZk requirements. Here, the and momentum range is required to be the 
same as for the signal. Since the momentum and angular distributions are slightly different 
for D° data and signal MC, the KID efficiencies are reweighted as a function of the polar 
angle of the signal track with respect to the e~ beam. In addition to the KID requirement, 
positively identified electrons are rejected using a similar likelihood ratio that also includes 
the energy deposited in the ECL. 

The dominant background is due to the e + e~ — > qq (q — u, d, s, c) continuum processes. 
A large MC sample shows that backgrounds from the b — > c transition are negligible since 
the momenta of their decay products are smaller than those in the signal decays. On the 
other hand, the momenta of the decay products from b — > u, s, d transitions other than 
the signal (denoted as other charmless B decays) can be as large as those in the signal 
decays. Events from these charmless B decays populate the negative AE region because of 
the energy carried away by a photon or 7T meson, which is not used in the B reconstruction. 
We take these events into account in the signal extraction as discussed later. 

We discriminate signal events from the qq background by the event topology. This is 
quantified by the Super- Fox- Wolfram (SFW) variable 0], which is a Fisher discriminant [l8| 
formed from modified Fox- Wolfram moments 19]. The angle of the S-meson flight direction 
with respect to the beam axis in the CM frame (9 b) provides additional discrimination. A 
signal likelihood ratio 1Z S = C S /(C S + £ q q) is used as the discriminating variable, where C s 
denotes the product of the individual SFW and 9b likelihoods for the signal, and C qq is that 
for the qq background. The probability density functions (PDFs) used for the likelihoods 
are derived from the MC for the signal, while events in the Mb c sideband (5.2 GeV/c 2 < 
Mb c < 5.26 GeV/c 2 in the AE acceptance) are used for the qq background. We make a 



mode-dependent requirement on 7Z S that maximizes N^ xp /^j iVf xp + N^ p , where iV| xp and 
A^| p denote the expected signal and qq yields based on our previous measurements [f| (upper 
limits are used for KK modes). The 1Z S requirements eliminate more than 90% of the qq 
background for the signal efficiencies given in Table |U 

Signal yields are extracted using a binned maximum-likelihood fit to the AE distributions 
after all the event selection requirements discussed above. The fitting function contains com- 
ponents for the signal, qq background, and other charmless B decays. If applicable, possible 
reflections due to the K ± /n ± misidentification are included as additional components. All 
of the fit parameters other than the normalizations are fixed. The signal PDFs are based 
on the MC. For the modes with a ir° meson, the PDF is modeled with an empirically de- 
termined parametrization For the other modes, the sum of two Gaussian distributions 
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with a common mean is used for the PDF. Due to the 7r ± mass assumption, each meson 
in the final state results in a shift in the peak position of about —45 MeV. Discrepancies 
between the peak positions in data and MC are calibrated using B + — > D°n + decays, where 
the D° — > fr + 7r~7T° sub-decay is used for the modes with 7r° mesons and the D° —>■ K + n~ 
sub-decay is used for the other modes. Here, the same analysis procedure used for the signal 
is applied except for the daughter particle reconstruction. The MC-based AE resolutions 
are calibrated using invariant mass resolutions of high-momentum inclusive D decays. We 
use D° -> K~ir + for the B° -> K + ir~, n+n', and K+K~ modes, D + -> K s tt + for the 
B + -»■ K s ir + , K° S K+ and B° -> K S K° S modes, and D° -> R-n+n for the modes with a tt° 
meson. The momentum ranges and reconstruction procedures for the D daughter particles 
are required to be the same as those for the signal daughter particles. The signal PDFs are 
also used for the reflections. The PDF of the qq background is determined from the M^ c 
sideband data and modeled with a first-order polynomial for the K + K~ and K S K S modes 
and a second-order polynomial for the other modes. The PDF for the other charmless B 
decays is taken from a smoothed histogram of a large MC sample. The K + tt° and 7r + 7r° 
modes are fitted simultaneously with a fixed reflection-to-signal ratio that is determined 
from the measured KID efficiencies and misidentification rates. For other modes, all the 
normalizations are floated. All fit results are shown in Fig. ^ 

The obtained signal yields are listed in Table ITT1 together with their statistical significances 



S = J— 2 ln(£ /AvJ, where C and Cn s denote the maximum likelihoods of the fits without 
and with the signal component, respectively. The fitted reflection yields are consistent 
within statistics with the expectations, which are derived from the fitted signal yields, KID 
efficiencies and misidentification rates, and the efficiencies of the R s requirements. We 
observe clear signals for B — > Ktt, 7r + 7r~ and 7r + 7r° decays. For the decays B — > KK, 
no significant signal is observed. We apply the Feldman- Cousins frequentist approach with 



systematic uncertainties taken into account 21] to obtain upper limits on the yields at 



the 90% confidence level (CL); these are used to set branching fraction upper limits. The 
branching fractions and upper limits are listed in Table |H] Here, our latest measurement 
for B° — > n° n° [l3| based on a data sample of 152 million BB pairs is also listed for 
completeness. The hierarchy of the branching fractions, B(B — > Ktt) > BiB — > titt), is 
confirmed. More statistics are needed in order to firmly establish the position of B — > KK 
in this hierarchy. 

The systematic errors in the branching fractions are the quadratic sums of the systematic 
errors in the signal yields, uncertainties in the reconstruction efficiencies, and the 0.6% 
error in the number of BB pairs. The systematic errors in the signal yields come from the 
uncertainties in the fit procedure. In order to study the sensitivity to the signal and qq 
background PDFs, each shape parameter is independently varied by its error in the fit. The 
sensitivity to the contribution from other charmless B decays is evaluated by changing the 
minimum AE requirement to —100 MeV (—150 MeV) for the modes without (with) a 7r° 
meson, to exclude most of these events from the fit. The resulting changes in the signal yield 
are added in quadrature and assigned to the systematic errors on the signal yields as listed 
in Table |H] The uncertainties in the reconstruction efficiencies are listed in Table IIHI along 
with the test samples that are used. The uncertainty for the track finding efficiency in the 
high-momentum region is obtained by comparing the ratio of yields of fully reconstructed 
and partially reconstructed test samples in data and MC. The uncertainties in the Kg and 
7T° reconstruction efficiencies are obtained from similar comparisons of yield ratios in test 
samples. Here, the test samples are restricted to the same K s and ir° momentum ranges as 
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the signal. The experimental errors in the branching fractions of these decays |22j are added 
in quadrature. The uncertainties in the KID efficiencies and misidentification rates are due 
to the statistics of the data test sample. We also checked the effect of the difference in the 
hadronic environment between the signal and test samples. No significant effect is seen in 
the efficiencies and misidentification rates. The TZ S requirement for each mode is applied to 
data and MC test samples, and the difference is included in the systematic error. 

To a good approximation in the Standard Model, the relative weak phase between the 
penguin and tree amplitudes in Kit modes is 03. It is in principle possible to extract 03 if 
the hadronic uncertainties are under control. Several approaches to constrain 03 have been 
proposed using the ratios of partial widths for Ktc and tttt modes with model-dependent 
assumptions on the hadronic uncertainties p; the ratios give cancellations of these uncer- 
tainties. We calculate such useful partial width ratios as listed in Table HVl Here, the ratio 
of charged to neutral B meson lifetimes t b +/t b o = 1.083 ± 0.017 [22] is used to convert 
the branching fraction ratios into partial width ratios if necessary, and the total errors are 
reduced because of the cancellation of the partially common systematic errors. Applying 
the approach of Buras and Fleischer [23j], for illustration, our T(K + tc~) /2T(K°7i°) measure- 
ment excludes the region 29° < 03 < 83° at the 90% CL based on MC pseudo-experiments 
while that of 2T(K + 7i )/T(K 'n- + ) gives no constraint. These results are obtained without 
any assumption on the tree-to-penguin amplitude ratio, but neglecting re-scattering effects 
and taking the size of the electroweak penguin as in Ref. j23[. Although a more aggressive 
constraint on 3 can be derived by introducing further model-dependent assumptions on 
the hadronic uncertainties, a coherent study of these approaches is required to reduce the 
model-dependence on hadronic uncertainties and to determine 03. 

A naive expectation for the tree-dominated 7t + tt~ and 7r + 7r° modes predicts 
2r(7r + 7T°)/T(7r + 7r~) = 1. The deviation of our result from this expectation, as given in Ta- 
ble IIV1 is consistent with our previous measurement and would indicate the existence of a 
significant penguin contribution in the 7r + 7r~ mode if the color-suppressed tree contribution 
plays a minor role |24| . This penguin contribution complicates the extraction of 02 from the 
time-dependent CP asymmetry in the -it + it~ mode (referred to as " pen guin pollution") (25^ . 
Applying an approach based on the isospin relations in tttt modes [261 ] . our measured ratios 
r(7r+7r°)/r(7r + 7r-) and r(7r°7r°)/r(7r + 7r^) in Table \N\ give the 90% CL bound on the size of 
the "penguin pollution" \6\ < 56°; the CL is derived from MC pseudo-experiments. Here, 
we also use the partial-rate CP asymmetry in the 7r + 7r~ mode, A nn = +0.58 ± 0.15 ± 0.07, 
which is reported in Ref. 0. 

In conclusion, we have measured or constrained the branching fractions for the B — > Kn, 
7r + 7r _ , 7T + 7r° and KK decays with 85.0 million BB pairs collected on the T(4S) resonance at 
the Belle experiment. We observe clear signals for B — > Kit, ir + 7T~ and tt + tt decays and set 
upper limits on B —>■ KK decays. The hierarchy of branching fractions reported in earlier 
measurements is confirmed. These results have significantly improved statistical precision 
compared to our previous measurements and supersede them. The results can be used to 
give model-dependent constraints on 3 , as well as limits on the hadronic uncertainty in the 
time-dependent analysis of 2 . 

We wish to thank the KEKB accelerator group for the excellent operation of the KEKB 
accelerator. We acknowledge support from the Ministry of Education, Culture, Sports, Sci- 
ence, and Technology of Japan and the Japan Society for the Promotion of Science; the 
Australian Research Council and the Australian Department of Education, Science and 
Training; the National Science Foundation of China under contract No. 10175071; the De- 



7 



partment of Science and Technology of India; the BK21 program of the Ministry of Education 
of Korea and the CHEP SRC program of the Korea Science and Engineering Foundation; 
the Polish State Committee for Scientific Research under contract No. 2P03B 01324; the 
Ministry of Science and Technology of the Russian Federation; the Ministry of Education, 
Science and Sport of the Republic of Slovenia; the National Science Council and the Ministry 
of Education of Taiwan; and the U.S. Department of Energy. 



[1] K. Abe et al. (Belle Collaboration), Phys. Rev. D 66, 071102(R) (2002). 
[2] B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 89, 201802 (2002). 
[3] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963); M. Kobayashi and T. Maskawa, Prog. Theor. 
Phys. 49, 652 (1973). 

[4] See, for example, the review "CP Violation in B Decay — Standard Model Predictions" by 
H. Quinn and A. I. Sanda in Ref. j^]. (a(= fa) and 7(= 03 ) are also used in the literature.) 
[5] See, for example, the chapter "Future Directions" and its references, in "The CKM Matrix 

and the Unitarity Triangle", M. Battaglia et al. (ed.), |hep-ph"7"0 304132 (2003). 
[6] B. C. K. Casey et al. (Belle Collaboration), Phys. Rev. D 66, 092002 (2002). 
[7] Y. Unno et al. (Belle Collaboration), Phys. Rev. D 68, 012001 (2003). 
[8] K. Abe et al. (Belle Collaboration), Phys. Rev. D 68, 012001 (2003). 
[9] K. Abe et al. (Belle Collaboration), |hep-ex/0401029| (2004). 
[10] B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 87, 151802 (2001); Phys. Rev. 

Lett. 91, 021801 (2003). 
[11] B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 89, 281802 (2002). 
[12] D. Cronin-Hennessy et al. (CLEO Collaboration), Phys. Rev. Lett. 85, 515 (2000); S. Chen 
et al. (CLEO Collaboration), Phys. Rev. Lett. 85, 525 (2000); D. M. Asner et al. (CLEO 
Collaboration), Phys. Rev. D 65, 031103(R) (2002). 
[13] S. H. Lee et al. (Belle Collaboration), Phys. Rev. Lett. 91, 261801 (2003). 
[14] B. Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 91, 241801 (2003). 
[15] A. Abashian et al. (Belle Collaboration), Nucl. Inst. Meth. A 479, 117 (2002). 
[16] S. Kurokawa and E. Kikutani et al., Nucl. Inst. Meth. A 499, 1 (2003) and other papers in 
this volume. 

[17] R. Brun et al, GEANT 3.21, CERN Report No. DD/EE/84-1 (1987). 

[18] R. A. Fisher, Annals of Eugenics 7, 179 (1936). 

[19] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[20] J. E. Gaiser et al, Phys. Rev. D 34, 711 (1986). 

[21] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998); J. Conrad et al, Phys. Rev. 

D 67, 012002 (2003). 
[22] K. Hagiwara et al. (Particle Data Group), Phys. Rev. D 66, 010001 (2002). 
[23] A. J. Buras and R. Fleischer, Eur. Phys. J. C 16, 97 (2000). 
[24] M. Gronau and J. L. Rosner, Phys. Rev. D 65, 013004 (2001). 
[25] M. Gronau and D. London, Phys. Rev. Lett. 65, 3381 (1990). 
[26] M. Gronau et al, Phys. Lett. B 514, 315 (2001). 



8 



TABLE I: Signal efficiencies for kinematic reconstruction (Rec), TZk requirements and TZ S require- 
ment along with the sub-decay branching fraction (/3 su t>) for K° — ► — > it + tt~~ and total signal 
efficiencies. 



Mode 


Rec 


TZ K 


TZ S 


£>sub 


Total 




0.731 


0.769 


0.672 





0.378 


K+7T 


0.461 


0.844 


0.501 




0.195 


K°7T+ 


0.571 


0.911 


0.560 


0.343 


0.100 


K°7T° 


0.314 




0.673 


0.343 


0.073 


7T + 7T~ 


0.756 


0.830 


0.560 




0.352 


1T+1T 


0.476 


0.911 


0.395 




0.172 


K+K- 


0.727 


0.713 


0.387 




0.201 


K+K° 


0.539 


0.844 


0.388 


0.343 


0.061 


K°K° 


0.447 




0.561 


0.235 


0.059 



TABLE II: Signal yields (N s ), statistical significance (<S), and branching fractions (B) for the 
B — > Kit, tt + tt~, 7r + 7r° and KK decays. The first and second errors are the statistical and 
systematic errors, respectively. For completeness, the tt°tt results from Ref. [ill] are also listed. 



Mode 


N s 




S[a] 


B [10- 6 ] 


K + it- 


595.9 1 ti t 


7.8 
7.7 


24.1 


18.5 ±1.0 ±0.7 


K+TT° 


198.9 ±21.5 t 


15.6 
4.8 


10.8 


12.0 ± 1.3 i H 


K°7T+ 


187.0 ± 16.3 ± 


1.5 
1.7 


16.4 


22.0 ±1.9 ±1.1 


K°7T 


72.6 ± 14.0 + 


4.9 
5.5 


5.8 


11.7 ±2.3 1 \ J 


7T + 7T~ 


132.7 i Hi ± 


2.7 
2.9 


8.5 


4.4 ±0.6 ±0.3 


7T+7T 


72.4 ± 17.4 + 


3.7 
3.4 


4.5 


5.0 ±1.2 ±0.5 


7r°7r° 


25.6 1 Ht 


1.6 
1.4 


3.5 


1.7 ±0.6 ±0.2 


K+K- 


-i-o t ti 




0.0 


< 0.7 


K+K° 


8.6 ± 5.9 




1.6 


< 3.3 


K°K° 


2.0 ± 1.9 




1.3 


< 1.5 



TABLE III: Uncertainties on the reconstruction efficiencies (<5e) along with the test samples that 
are used. 

Source 8e/e [%] Test sample 

track finding 1.0 D*+ -> D°(^ K%ir + ir-)ir + 

K% 4.4 D + -» K%n+, K~Ti + T: + 

vr° 3.5 7] -> 7r 7r°vr , 77 

TZ K 0.2 D*+ -> D°(-> K--k+)tt + 

TZ S 1.3-7.8 B + -> D°(-> K+7T-, K+tt-tt )^ 
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TABLE IV: Partial width ratios of B — > Kir and tttt decays. The errors are quoted in the same 
manner as in Table ILTI 

Modes Ratio 

T(K + ir-)/ T(K°7r + ) 0.91 ± 0.09 ± 0.06 

T{K + tt-)/2T(K°tt°) 0.79 ± 0.16 ± 0.09 

2T(K + tt )/ T{K tt + ) 1.09 ± 0.15 ± 

r(vr + vr-) / T(K+ir-) 0.24 ± 0.03 ± 0.02 

r(vr + vr ) / T(K°tt°) 0.39 ± 0.12 ± 0.06 

2r(vr + vr ) / r(K°ir+) 0.45 ± 0.12 ± 0.05 

2r(vr + 7r°) / T(tt + tt-) 2.10 ± 0.58 ± 0.25 

r(vr 7r ) / r(vr+7r-) 0.39 ± 0.15 ± 0.05 

r(vr 7r ) / r(vr + 7r°) 0.37 ± 0.16 ± 0.05 
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FIG. 1: AE distributions for B — ► Kir, 7t + n~, 7r + 7r° and KK decays. Fit results are shown 
as the solid, dashed, dotted and dot-dashed curves for the total, signal, qq background and the 
other charmless B decays, respectively. In addition, reflections due to K^/ir^ misidentification 
are shown as hatched areas. 
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